Symmetry constraints for dispersionless integrable 
equations and systems of hydrodynamic type 

L.V. Bogdanov* and B.G. Konopelchenko^ 



Abstract 

Symmetry constraints for (2+l)-dimensfonal dispersfonless inte- 
grable equations are considered. It is demonstrated that they naturaUy 
lead to systems of hydrodynamic type which arise within the reduction 
method. One also easily obtaines an associated complex curve (Sato 
function) and corresponding generating equations. Dispersionless KP 
and 2DTL hierarchy are considered as illustrative examples. 



1 Introduction 

Dispersionless integrable equations arise in various contexts and have va- 
riety of applications from complex analysis to topological field theory (see 
[T]-j22j). The reduction method seems to be the most developed method of 
construction of exact solutions of dispersionless equations [Bj l6l Isj [T7].[23]- 
It reveals a close connection between such equations and certain sys- 
tems of hydrodynamic type, emphasising a role of complex curves (Sato 
functions). The reduction method has been successfully applied to several 
dispersionless equations. However, an intrinsic background of its validity 
and efficiency seems to remain not very clear. 

In the present paper it is shown that systems of hydrodynamic type 
arising within the reduction method are nothing but the symmetry con- 
strained (2-|-l)-dimensional integrable equations. We will consider here the 
dispersionless KP and 2DTL hierarchies. 

First, we find a general infinitesimal 'nonisospectral' symmetry transfor- 
mation of the dispersionless r-function (free energy F). Then we impose 
symmetry constraint of the type 'isospectral symmetry = nonisospectral 
symmetry', similar to dispersionfull case |35j-j4Uj. It is demonstrated that 
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such a constraint reduces (2+l)-dimensional dispersionless equation to a set 
of (l+l)-dimensional systems of hydrodynamic type with a finite number 
of dependent variables. We present also a simple way to find an associated 
complex curve and to construct generating equations. 

Our approach starts with the (9-method applied to dispersionless equa- 
tions. It provides us with a simple and natural way to consider symmetries 
and symmetry constraints. At the same time, an analysis of such constraints 
indicates a direction of a possible extension of the quasiclassical 9-method. 

2 Quasiclassical 5-method 

Dispersionless integrable hierarchies can be described in different forms 
within different approaches. In the papers j41l 1421 144j it was shown that 
such hierarchies can be introduced starting with the nonlinear Beltrami 
equation (quasi-classical 9-problem) 

S, = W{z,z,S,), (1) 

where z € C, bar means complex conjugation, Sz = Sz = ^^gj^j 

and W (quasi-classical 9-data) is an analytic function of Sz, 

oo 

W{z,z,Sz) = Y.Wp{z,z){SzY. (2) 

p=0 

Applying the quasi-classical 9-dressing method based on equation (P), one 
can get dispersionless integrable hierarchies and the corresponding addition 
formulae in a very regular and simple way. Such an approach reveals also the 
connection of dispersionless hierarchies with the quasi-conformal mappings 
on the plane. 

The specific choice of the hierarchy depends on the choice of domain G 
which is the support of the (?-data. Once the domain is fixed, we introduce 
the class of functions Sq{z), having the derivative dzSQ{z) analytic in G, and 
formulate the boundary problem in G in the following way. Let the function 
Sq{z) be given. The problem is to find the function S = Sq + S , satisfying 
(jTj) in G, with S analytic outside G and decreasing at infinity 

Introducing parameterization of the class of functions 5*0(2) in terms 
of infinite number of variables (times), and using the technique of quasi- 
classical 9-dressing method, it is possible to demonstrate that S{z,z,t) is a 
solution of Hamilton-Jacobi equations of corresponding dispersionless hier- 
archy. 
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For dispersionless KP (dKP) hierarchy we take unit disc D as G, and 
parameterization of 5*0(2:) is given by 

00 

5o(z,t) = ^t„z", (3) 

n=l 

S{z, z, t) = Sq + S , S is analytic outside the unit disc, and at infinity it 
has an expansion S = Yl^i Si{t)z~'^ . The quantity p = ^ is a basic 
homeomorphism j42j . Important role in the theory of dKP hierarchy is 
played by the equation 

p(z)-p(zi) + ziexp(-D(zi)5(z)) =0, zGC, 2iEC\Z), (4) 

(where D{z) is the quasiclassical vertex operator, D{z) = Yll^=i'k,^^^ 
\z\ > 1) which generates Hamilton- Jacobi equations of the hierarchy by 
expansion into the powers of z^^ at infinity (see, e.g., |47j ) . This equation 
also implies existence of the r-function, characterized by the relation 

S{z,t) = -D{z)F{t), (5) 

and provides the dispersionless addition formula. 

For the dispersionless 2DTL (d2DTL) hierarchy the (9-data are localized 
on the domain G which is an annulus a < |z[ < b, where a,b {a,b G M, 
a,b > 0; b > a) are arbitrary, Sq can be represented as jHl El HZj 



00 



Soit, X, y) = t log Z + ^ Z^Xn + ^""Vn- (6) 

n=l n=l 

We assume that S{z) ~ Yl'^=i as z ^ oo and denote 5(0) = (p, G+ = 
{z, \z\ > b}, G- = {z, \z\ < a}. The functions p+ = and p_ = ^ 
have pole singularities while p = ^ has a logarithmic singularity. Relations, 
characterizing the r-function F, are (p = DF, S{zi) = —D+{zi)F (z\ e G+), 
5(Z2) = 4>- D^{Z2)F (Z2 e G_), where D+{z) = ^Zi ll^i-^ D^{z) = 
X^nLi ^-^"a^; ^ ~ Hamilton- Jacobi equations of the hierarchy are 
generated by the relations |S1 ^1 

^P{z) _ gp(^i) ^ ^_^^-D+{z,)S{z) z G C, zi G G+, (7) 

e^'(^) (l - e-^-(-2)S(.)^ ^ z G C,Z2 G G_, (8) 

these relations also provide dispersionless addition formulae for F. 
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3 r-function and symmetries 

It was noted in 01j that equation (pQ) is a Lagrangian one. It can be obtained 
by variation of the action (for the boundary problem in G) 

f = -^Jj^ Q^J. - W.iiz, z, S,)^ dz A dz, (9) 

where 

W-i{z,z,S^) =y2wp{z,z) \ ^ , dr,W-i{z,z,r]) = W{z,z,7]). 

p=0 ^ 

One should consider independent variations of S", possessing required ana- 
lytic properties (analytic outside G, decreasing at infinity), keeping 5o fixed. 

In |48j it was proved for dispersionless KP and 2DTL hierarchies that 
the action evaluated on the solution of nonlinear Beltrami equation 
gives a r-function of the corresponding hierarchy. Here we formulate first 
an auxiliary statement, valid for an hierarchy with arbitrary G, and then 
recall the results of [15] . 

Lemma 1 Let us evaluate the action 0) on the solution of the boundary 
problem for Q) and consider it as a functional of Sq, 

^^^^^ " "2^ IJg ~ ^' ^'^) '^^ ^ 

Then 

5F = -^ [ Sd6So. (11) 
2vri Jqg 

Using the formula (jllj) . it is easy to reproduce the statements of the work 

Proposition 1 The function 

F{t) = JJ^ Q5,(t)5,(t) - W^iiz, z, S,{t))^ dz A dz, (12) 

i.e., the action (0) evaluated on the solution of the problem ^\|, with Sq 
given by is a r-function of dKP hierarchy. 
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Proposition 2 The function 

F(t,x,y) = -^ // (ls2{t,x,y)S,it,x,y) -W.i{z,z,S,{t,x,y))] dz Adz, 



27ri JJg \2 

with Sq given by evaluated on the solution of the problem. 0), is a r- 
function of dispersionless 2DTL hierarchy (G is an annulus defined above). 

The function W is the d data for the quasiclassical 5-problem. Its variations 
provide us with a wide class of variations of the function F (first we will use 
formula defined for arbitrary G). 

For the functions W of the form @, varying vun{z, z), one has 

oo oo i- 

6w = Y, 5y^n{z, -z) {s.T , sw., = Y.^ {s.r^' , 



.n+1 

n=l n=l 



and 



27ri 



SF = ^ II U{z,z,S^)dz Adz, (13) 

IG 



where U {z, z, r\) is an arbitrary function with the support in z-plane be- 
longing to G, analytic in r]. The formula (|13j) gives a general variation 
(infinitesimal symmetry transformation) of the r-function. Let us consider 
some special symmetry transformations. 

Considering elementary variation ^Wn^ = eanoS{z — zq), Swn = 0,n ^ uq, 
one gets 



2m (no + 1) ° 



and, respectively, 



Taking superposition of elementary variations H14|) , we obtain a variation of 
the form 

6F = ef{S.{zo)), (16) 

where / is an arbitrary analytic function (summation over different points 
and integration over zq are also possible). 

Another simple symmetry transformation is given by 

6F = ec{S{z2) - S{zi)), (17) 
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where zi, Z2 are arbitrary points belonging to the domain G, and c is some 
constant. Indeed, 



thus the function S{z2) — S{zi) corresponds to the special case of general 
variation In the limit Z2 — > zi the transformation H17|) is reduced to 

elementary symmetries dJ). 

4 Symmetry constraints 

Given infinitesimal symmetries of the form 

6F = e$, 

(cf., (fn^ . ifTBl) . ((T7|) l. it is possible to introduce symmetry constraints 



These constraints are preserved by the flows of the hierarchy and represent 
a reduction of the hierarchy. There is a lot of similarity between constraints 
H18() and symmetry constraints in the (2+l)-dimensional dispersionfull case, 
which are rather well studied |35j-|4f)j. Below we will demonstrate that 
constraints (jlHj) lead to hydrodynamic reductions, and give a regular way to 
introduce basic objects connected with these reductions. 

Technically, in the dispersionful case one starts with linear equations 
of the hierarchy (which give the hierarchy itself in the form of compati- 
bility conditions). Symmetry constraints connect the potentials in these 
equations with some special wave functions, and in the end one obtains 
(l-l-l)-dimensional integrable system for the set of wave functions. The sim- 
plest symmetry constraint for KP hierarchy leads to nonlinear Schrodinger 
equation. 

In the dispersionless case we will follow the same strategy. The role of 
linear equations of the hierarchy is played by Hamilton-Jacobi equations of 
the dispersionless hierarchy. 





(18) 
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5 Zakharov reduction for the dKP hierarchy 



We begin with a simple example of symmetry constraint for the dKP hier- 
archy. The first Hamilton-Jacobi equation of dKP hierarchy (dispersionless 
analogue of Lax equation) is 

Sy=p^ + 2u; p = S,,u = -dJi, (19) 

where x = ti, y = t2. Let us consider the symmetry constraint 

F,. = f{S,{zo)), (20) 

which is equivalent to the condition (cf., ©) 

u = dxf{Sz{zQ)). 

Differentiating ()19() with respect to z and evaluating the result aX, z = zq, 
we get 

/c, / ^ = 2po, Po = P{zo). 

<Jx^z{Z0) 

Using this equation and evaluating (|19jl at z = zqi we obtain a system of 
hydrodynamic type 

r dyPo = d^{pl + 2u), 
\ dyU = 2d^{pou). 

This system is well known |2|, it represents a dispersionless limit of nonlinear 
Schrodinger eqution. Thus a constraint (|2()|) represents a proper dispersion- 
less analogue of the simplest symmetry constraint of KP hierarchy. 

It is possible to move further and demonstrate that the symmetry con- 
straint H2U|) leads to explicit expression for the function z{p), which plays 
a fundamental role in the picture of constrained hierarchy and allows to 
find Riemann invariants and hodograph equations. To do this, we will use 
generating Hamilton-Jacobi equation (jlj) rewritten as 

D{z)S{z') = - log P(^)-P(^'\ > 1^ (22) 

z 

which produces infinite set of Hamilton-Jacobi equations by expansion into 
the powers of z~^ at infinity (see, e.g., |3Z1)- Differentiating this equation 
by z' and evaluating the result at z' = zq, we get 

D{z)Sz{zo) ^ 1 
dxSz{zo) p-po' 
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Then, using the constraint ()2U() and simple formula 

D{z)F^ = z-p, 

we obtain the relation 

z = p + ^^. (23) 
P-Po 

The function z{p) (Sato function) plays a crucial role in the picture of dis- 
persionless hierarchy. Dispersionless KP hierarchy (for general z{p)) can 
be written in the form dnz{p) = {z",^}. Relation ((231) defines the function 
z{p) for the constrained hierarchy, it is well known 2 and is usually called 
Zakharov reduction. Thus we have demonstrated that symmetry constraint 
(|20j) leads to Zakharov reduction (|23j) . 

Riemann invariants for the constrained hierarchy of equations of hydro- 
dynamic type (the first of which is pT]) ). are given by the values of z{p) ((23|) 
at Pi where ^ vanishes, Xi = z{pi), f|(Pi) = 0. It is rather straightforward 
to introduce hodograph transform for the solution of constrained hierarchy 

6 Generic symmetry constraint for the dKP hier- 
archy 

Now we will use the symmetry (|17() and consider the constraint 

N 

F, = Y,c^{Si-S,), (24) 

i=l 

where Si = S{zi), Si = S{zi), Zi, Zi are some sets of points, and Cj are 
arbitrary constants. In equivalent form, 

N N 

u = dx'^Ci{Si - Si) = '^Ci{pi -pi). 

i=l i=l 

Evaluating Hamilton-Jacobi equation ((T^ and its higher counterpart 

St=p^ + Zup-^Si, 

where t = ta, at z equal to Zj, Zj, we immediately get two systems of hydro- 
dynamic type 

dyPk = dx{(pl) + 2 ^ Ci(pi - Pi)) 

^ k = l,...,N (25) 

dyPk = dx [{Pk) + 2 2_^Ci{pi- Pi)) 
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and 

dtPk = dx{{pi) + 3pk Ci{pi - Pi) + 2 ~ )) 

i i 

dtPk = dx {{Pk) + 3pk Ci{pi - Pi) + 2 X] ~ Pi^) 



(26) 



These systems can be written in Hamiltonian form, 



H3 
^ 3 



H4. = \j dx Ci(p^ - Pi^) + 6 ^ Ci{pi - Pi) X ^ Ci(p2 - j ^ 

where p = (pi, . . . ,pn,Pi, ■ ■ ■ iPnY, C is diagonal N x N matrix with the 
entries q, Cik = CiSik, i,k = 1, . . . ,N. 

Common solution of the systems (|25|). (|2H|) generates a solution u = 
J2i (^i{Pi ~ Pi) of the dKP equation. 

The Sato function z(p), which allows to construct Riemann invariants 
and higher Hamiltonians, is the central object for the constrained hierarchy. 
To find it, we, similar to the first example, use the generating Hamilton- 
Jacobi equation (j^H). Evaluating this equation at z' equal to Zi, Zi and 
combining the results, we get 

AT AT 

D{z) V Ci{S, -~S,) = -y^c, log (27) 



Then, due to the constraint H24|l . 

D(z)F,. = -J^Cilog 

i=l 

and finally we obtain the function z{p), 



N 

p-Pi 



,=1 V-P^ 



z = p-J2c^log^. (28) 

U P'P' 

Equation |^ = is an algebraic one, having 2N roots pa, and Riemann 
invariants for the constrained hierarchy are Aq, = z{pa). 
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Expansion of z(j)) at p ^ oo looks like 



oo N 



P 



n=l ^ 1=1 



Equation ()22() allows us to obtain the generating equation for the whole 
hierarchy of systems of hydro dynamic type associated with the constraint 
(|24() . Indeed, differentiating H22|) with respect to x and evaluating it at the 
points Zi, Zi, one gets 

Diz)p, = -dMp-P>.) , = (30) 

where p is a function of z, {pi, . . . ,p]y,pi, . . . ,p]y), defined by the relation 
(|28|) . Expanding both sides of this system into the powers of z~^, one gets 
the systems (|25|) . (j26p and their higher counterparts. Introducing to So{z, t) 
dependence on extra time tiog, So{z') SQ{z') — tioglog{l — gf— = D{z), 
we may consider pOf) as a system of hydrodynamic type with the time tiog- 
Hamiltonian for this system is given by the expression 



Hiog = / -p^ - V Q {pi -Pi+Pi log(p - Pi) - Pi log{p - Pi)) . 




Differentiating H30|) by z, we will get another form of generating system, 

Dz{z)pk = -dx 
Dz{z)pk = -Ox — 

P-Pk 

where p,pz are functions of z, (pi, . . . ,pN,pi, ■ ■ ■ ,Pn), defined by the relation 
H28() . This system can be considered as a system of hydrodynamic type with 
the time tpoio So{z') So{z') + tpoieiz — z')~^ , ^ = Dz{z). Hamiltonian 
for the generating system H31() is remarkably simple, 

^^poie = j pdx, 

and its expansion into the powers of z~^ gives a general formula for the 
Hamiltonians of constrained hierarchy, 

Hn = ^j reSoo(^(p)") dx = dJn, (32) 
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which is in agreement with general results 

In a similar way one can treat the symmetry constraints 



= ^CiiS.-S.). (33) 

1=1 



In particular, since 



1^ = + n„_2P" 2 + • • • + no, 
otn 

one obtains a corresponding Sato function 

AT _ 

E{p) =z-=p- + n„_2p"-2 + • • • + no - V Q log . (34) 

Remark. It is easy to check that in the limit Zi = Zi + e.; — > 0, 0^6.; = const, the 
Sato function reproduces the curve 



N 

2"=p" + n„_2P"-' + --- + no + y (35) 

which has been discussed in ^1 El ■ At n = 1 one gets a general Zakharov reduc- 
tion. 

7 Symmetry constraint for the d2DTL hierarchy 

Now we will consider symmetry reduction of d2DTL hierarchy defined by 
the constraint 

N 

Ft = Y,Ci{Si-Si). (36) 

i=l 

Technically this case is very similar to that considered in the previous sec- 
tion, and we will omit some details. The first Hamilton-Jacobi equations of 
d2DTL hierarchy are 

d,S{z) = eP - dtSi, (37) 
dyS = e**-P, (38) 

where x = x\, y = yi. Using the constraint H36|) . from these equations we 
get two systems of hydrodynamic type, namely 

<^ ^ ; - ^ ^ - ( i = l,...,iV 39 
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and 

( dyPi = dtexp{-pi + J2kCk{Pk-Pk)) .^^ ^ ^^Q^ 

\ dyPi = dtexp{-pi + Y^f,CkiPk-Pk)) 

Both systems can be written in Hamiltonian form, 

,.p..^, a,,./^, ^-{^^{V -c-))*- 

where E is 2N x 2N matrix with all entries equal to 1, Eik = 1, i,k = 
1,...,2N, C is diagonal N x N matrix with the entries q, Cik = Ci5ik, 
i,k = 1, . . . , N , and the Hamiltonians are, respectively. 



^1= dtJ2ck{eP'-e^')^ (41) 
k 

H{ = Jdt exp CkiPk - Pk?j ^ Ckie-P" - e-P>'). (42) 

Using generating equation 0, we obtain the function z{p) for the con- 
strained hierarchy. 



k=l 



n(^) ■ («) 



Starting with second generating equation (jH)), we get equivalent form of this 
relation, 

z = exp{p-}_^Ck{pk-Pk)) 11 (^ g-p_g-p, J ■ (44) 

Equation |^ = is algebraic (with respect to the variable ^ = e^), having 
2N roots and Riemann invariants for the constrained hierarchy are = 

Generating systems for the constrained hierarchy read 



D+iz)pk = -dtlog{eP - eP") 
D+{z)pk = -dtlog{eP - eP'^) 

D^{z)pk = -dtlog{l - eP-P'^] 
D^{z)pk = -dtlog{l - eP-P'^] 



k = l,...,N, (45) 
k = l,...,N, (46) 
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('logarithmic' form) and 

(Dz±{z)pk 
Dz±{z)pk 



ePk-p - 1 

Pz 



ePk-p - 1 



Pz 



k = l,...,N, 



(47) 



('pole' form). Hamiltonian for the 'pole' form of generating equations i 

H±{z) = - I pdt. 



Expansion of this Hamiltonian at z = oo and z = gives general expressions 
for the Hamiltonians of constrained hierarchy, 



where ^ = e^. 

In more detail, the Hamiltonian and bi-Hamiltonian structure of equa- 
tion considered above will be considered elsewhere. 

Symmetry constraints of the type Fx„ = '^f=iCi{Si — Si) and Fy^ = 
J2iLi^ii^i — Si) are treated analogously. The corresponding Sato function 
is given by (|5H) with substitutions p ^ e^, p ^ e~^, respectively. 

8 Conclusion 

As we have seen, the systems of hydrodynamic type and associated complex 
curves (Sato functions) are direct consequences of the generating Hamilton- 
Jacobi equations (|7I8() under the symmetry constraint. 

We would like to note here that complex curves discussed above are con- 
nected with the 9-equation in a very natural way. Let us consider the dKP 
hierarchy as an example. The function p = dxS solves a linear Beltrami 
equation pz = W'pz, and under certain conditions is a basic homeomor- 
phism. The inverse function z{p,p) is analytic in some neighborhood of 
infinity. The Cauchy- Green formula implies that 





(48) 



(49) 




(50) 



In the neighborhood of infinity one has 



oo 




(51) 
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where 

1 ff rn-l dz 



2niJJ^'"'W'^^'^'^^'' ^^^^ 

Let us choose a special || with a support along the curve F = |J • Fj com- 
posed of disconnected pieces Fj with the ends at the points Zi and Zi and 
such that || = — X^i^iCjJrj- For such 'constrained' 9-data formula (jSU]) 
gives 

N 

P - 2^ Ci log 



z 



while from (|52() we obtain 



n 

that coincides with (|28I29|) . For the 3-data given by || = — "^^^i ai6{p—pi), 
equation (|5fl|) immediately leads to the Sato function associated with Za- 
kharov reduction (see (|H5|) ). Thus, the Cauchy-Green formula is a generator 
of complex curves associated with equations of hydrodynamic type. 

Equation is the two-dimensional extension of the integral equation 
discussed in [Ij, |^. 
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